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ABSTRACT 



We study the phase structure of full QCD within the canonical ensemble 
(CE) with respect to triality in a lattice formulation. The procedure for the 
calculation of the effective potentials (EP) in the CE is given. We calculate 
the EP for the three dimensional SU{2) gauge model at finite temperatures in 
the strong coupling region. The potential exhibits a genuine deconfinement 
phase transition unlike the similar potential obtained in the grand canonical 
ensemble (GCE) which demonstrates explicit Z{Nc) symmetry breaking at 
any temperature. Furthermore, we investigate the EP with the chiral con- 
densate included. In contradiction to other authors, we find chiral symmetry 
restoration in all triality sectors. In the scheme with massless staggered 
fermions we observe chiral symmetry restoration accompanying a deconfine- 
ment phase transition of first order. Above the critical point, besides two 
Z{2) symmetric " deconfining" vacua there exists a metastable "confining" 
vacuum in a wide region of the (Aj, 7)-plane. Such a picture could be in- 
terpreted as an indication on a mixed state of hadrons and quarks in the 
vicinity of the critical line. 
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A comprehensive understanding of the phase structure of Quantum Chro- 
modynamics (QCD) is the basis for many phenomenological apphcations, for 
possible signals of the quark-gluon plasma and other topics (see, for instance, 
the review on the Symposium 'Lattice-94' In numerical calculations a 
number of new results has been obtained. Nevertheless, the theoretical un- 
derstanding of the phase structure of full QCD, despite huge efforts, is almost 
on the same level as it has been around fifteen years ago. In our view, one 
of the central points of the problem is a reliable theoretical determination 
of the deconfinement phase transition in QCD with dynamical fermions. It 
was understood in the beginning of eighties that in pure gluodynamics 
deconfinement can be described in terms of the spontaneous breakdown of 
the global Z{Nc) symmetry of the vacuum. An appropriate order parameter, 
the Polyakov loop, was proposed to distinguish the confined phase of static 
colour charges from the deconfinement. Since loops of dynamical quarks vi- 
olate in general the Z{Nc) symmetry explicitly, a naive application of this 
picture to full QCD is not possible. At large values of quark masses a reminis- 
cent of deconfinement can be seen in the behaviour of some thermodynamical 
functions. In QCD with massless quarks the behaviour of the quark conden- 
sate demonstrates a phase transition associated with the restoration of chiral 
symmetry. It was stressed, however, that "with two flavors of light quarks in 
the staggered fermion scheme, there appears to be a dramatic crossover, but 
so far no evidence for a genuine phase transition" [|l| (see also for review 
of the thermal transition in QCD). 

Recently, in several articles 0, ^, |^ strong arguments have been given that 
the general statement "dynamical quarks in QCD breaks Z{Nc) symmetry 
explicitly", is not valid. It has been discussed that the fermionic determinant 
in finite temperature QCD describes a grand canonical ensemble (GCE). One 
can restrict the statistical averaging procedure to a canonical ensemble (CE) 
with n-ality zero which may still be the grand canonical with respect to 
baryon number. This CE does not break Z{Nc) symmetry explicitly. 

In Ref. P] we speculated about the possible QCD phase structure. We 
gave some arguments that there should be a genuine deconfinement phase 
transition while in the GCE only a crossover has been observed. Here, we 
would like to supplement our intuitive arguments by a calculation of the 
effective potential in the CE at strong coupling. Since the Z{N^ symmetry 
is not explicitly broken in this ensemble by dynamical quarks, one may hope 
to find a phase transition similarly to pure gluodynamics. Furthermore, we 
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calculate the effective potential with the chiral condensate included. We shall 
show that the chiral symmetry is restored in all n-ality sectors. This is in 
disagreement with the GCE description of full QCD Q and is another 
argument in favour of our approach. First, we review our general strategy 
and then apply it to (2+l)-dimensional QCD with Wilson fermions. Further, 
we use our approach to study the behaviour of chiral symmetry within the 
CE in the scheme with staggered fermions. 

Let us consider the effective model of the Polyakov loops (PL) W^g in the 
strong coupling region calculated first in Ref.0: 

where So is the pure gluonic action and 7 ~ (p')^'- ^'^^ ^ system with 
heavy quarks one can apply the hopping parameter expansion of the fermionic 
determinant in the scheme with Wilson fermions 

Zg = Det[D,y] ^ e'^Ss^f^^^+oe^^), (2) 

with the hopping parameter h = (2^^)^*- Combining (|l|) and (Q) one gets 
for the effective theory of PLs in the GCE with respect to triality 

Zgce = J ]JdW^exp[^J2 ^P^-^P^^+n + hJ2 SpW^]. (3) 

X x,n X 

It is well known that the model (|^) does not demonstrate a critical behaviour 
at any 7 and h. This is a consequence of the explicit violation of the global 
Z[N) symmetry in the GCE. 

Generalizing a simple technics for describing the system of zero triality 
[^, ^], one obtains the corresponding partition function in the CE 

1 f 

ZcE = ^ XI y n ^-^^ b LccLcS+n + hY^ ZkLg] , (4) 
k=l X x,n X 

where Lg = SpWg is the trace of the Polyakov loop matrix and Zk = exp[^k] 
fixes the phase sector. The integration in has to be done with the invariant 
measure for the SU (N) group. We discuss now the SU (2) gauge group. A 
naive effective potential for the PL in the CE for = 2 is of the form 

Veff = 2d-fL^ + ln(l -L^) + ^ ln(cosh NhL), (5) 
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where N is the number of lattice sites at a given time-shce. In the thermody- 
namical hmit we have the term | Lh \ in (^) which breaks Z{2) symmetry at 
any temperature. It follows that < L > always differs from zero in this ap- 
proximation. One could interpret this breaking as dynamical caused by the 
quark sea at low temperatures and as spontaneous at high temperatures 
Thus, if this is the case we have two different mechanisms of breaking. To see 
the critical behaviour we would need an order parameter (OP) which could 
distinguish between these mechanisms. The PL obviously is not suitable for 
this goal. 

The starting point of our calculations was the simple observation that at 
least in the zero temperature case (or in the limit of infinite coupling) the 
naive mean-field does not lead to the expected result: Z{Nc) should not be 
broken and CE and GCE descriptions have to coincide since in fact we have 
no ensemble at zero temperature but only a system in its ground state: The 
expectation value of the PL at zero temperature equals zero, whereas the 
naive effective potential as well as the GCE predicts Z{Nc) is broken and the 
PL differs from zero in this limit. 

In this respect the essential point is: 

In different phase sectors in (^) the minimum of the effective action is 
reached for different Z{Nc) configurations. This is very important in the 
disordered confinement phase where the PL may take any value. When we 
consider the different phase sectors on the same configurations it follows that 
only one of them will survive in the thermodynamical limit. However, if these 
configurations are indeed disordered, averaging over them should give the 
opposite result: reaching equal minima but for different configurations, two 
sectors contribute to the thermodynamical limit in the confinement phase, 
cancelling all noninvariant contributions. 

Hence, the main idea is to execute the summation over Z{N) configura- 
tions contained in the PL and only then to apply the mean field approxima- 
tion. Doing in this way we can reproduce the true minima of the effective 
potential in each triality sector. To accomplish this one uses the formula 

I dUG{U) = ^T.j dUG{zU), z G Z(iV,), (6) 

which is valid for any compact group and its central subgroup Z{Nc). This 
formula gives us a chance to verify whether the final result depends explicitly 
on the Z{Nc) configurations contained in the initial action Sq or not. If yes. 
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we have the Z{Nc) symmetry broken at all temperatures. In the other case, 
we can find a spontaneous breaking only at high temperatures. 

Thus, we start from the following general expression for the EP applying 
Eq.® 

1 r 



k=l{s^},s^GZ{N,)- 



SxL^ - NL 



(7) 



where we the contribution to V^ff at zero temperature. We choose the 
gauge f/o(x, t) = 1 for t G {1, 2, ..,Nt — l} and Uq{x, Nt) = Wx with diagonal 
Wx. 

The effective action Sef / in (|^ is defined as an integral over all space com- 
ponents of gauge fields with the Wilson fermionic determinant included. Its 
calculation is clearly far from our present ability. Using the strong coupling 
and the hopping parameter expansions as well as the mean field approxima- 
tion one arrives to the following effective potential in the CE 

1 Nc 

exp [7L^ -ReSa;S*+„ + hL^ ZkRes^ + Vim] , (8) 

x,n X 

with 7 and h defined in (0) and (|), respectively. Vjm stays for the invariant 
measure for the SU{N) group. The summation over the Nc n-ality sectors 
guarantees that only states contribute to the effective potential whose number 
of spins Sx is a multiple of N^. 

Eqs.(|^)-(P) are still valid in arbitrary dimensions and for any SU{N) 
group. To make some further progress we consider in the following the SU{2) 
gauge theory in d = 2 + 1 dimensions. Then, our effective model (P) is two 
dimensional. We have to solve the 2 — d Ising model with an effective coupling 
in an external field and to investigate the effective potential obtained. 

After some calculations we find the following bound for the EP 

V;5f(L) < Vim{L) + j^\nZj{^L) +lncoshLh + 

-^ln^5]cosh[ArtanhL/i(< s(0)s(r) >y/^]. (9) 
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Zi in (|^) is the partition function of the 2 — d Ising model without external 
magnetic field calculated at the effective coupling 7^ = ^l?' . 

We used the following inequality to calculate the sum over pair correlation 
functions^ 

2k 

r = E < n ^(^^) > < E < Ko)s(r) . (10) 

The EP in exhibits a very essential feature of the CE description, namely 
all odd correlation functions are identically cancelled due to the global sum- 
mation over the center elements Zk- Just these contributions lead to the 
explicit violation of the Z{2) symmetry in the GCE description. It should 
be stressed that the cancellation happens exactly and is not a consequence 
of the approximations used. 

To make the following arguments as transparent as possible and to pro- 
mote the analytical evaluations to the very end we use simple asymptotics 
for the 2 — d Ising model with an external field at small and large values of 
the effective coupling 7/, |T^. This leads up to a constant independent of L 
to asymptotic expressions for the EP of the form 

V^^ff{L,-fL < 1) = \n{l- L^) +\n cosh Lh + 2-fL'^{tanh LhY + 

'y^L^[l + 6(tanhL/i)2 - 7(tanhL/i)4] + 0(7=^) (11) 



and 



K^f (L, 7L » 1) = ln(l - L') + F{h). (12) 



The function F{h) is defined as 

F{h) = 27L2 + Xh + u^ exp(-2X/i) + 2u^ exp(-4X/i) + 
M^[-exp(-4X/i) +6exp(-6X/i) +exp(-8X/i)] +0(m^), (13) 

where u = exp(— 47/^^) and X =\ L \. 

Analyzing the EP ( [TI]) one concludes that at sufficiently small 7 and 
h the maximum of the EP is located at L = 0. As 7 and/or h increases 
two maxima with L = ± | L^ax 1 7^ appear dominating the EP. In this 
approximation it is easy to find a critical line in the (7 — h)-p\a.ne. We omit 

^In our terms the effective potential presented in |q] would correspond to F < N^'^ . It 



is obvious, however, that the bound (fO) is lower. 



6 



this simple calculation as further we study in details the more interesting 
case with the quark condensate included. The EP for 7^ 3> 1 (|12|) is also 
formally symmetric under Z{2) transformations as it should be for the CE at 
any temperature. However, because of the presence of the term Xh in (|1^) 
the EP develops two maxima at L 7^ for any h. Thus, Z{2) symmetry can 
be spontaneously broken at high temperatures indicating the deconfinement 
phase transition. Since we have two degenerate maxima for 7^1, domain 
walls might separate regions with positive and negative values of L, similar to 



pure gluodynamics. Concerning this question see, however, [|TT|. It is worth 
to mention that the form of the EP at large 7^ is close both to the EP in 
the GCE and to the EP in the CE calculated in . This means that the EP 
from 1^ is only valid in the deconfinement phase. 

Now, we would like to construct and to study an effective potential for 
Polyakov loops and for the quark condensate a in the CE. It has been argued 
in 1^ that the chiral symmetry is restored only in one n-ality sector (e.g., in 
the real phase for SU{3)) while it is still broken in other sectors above the 
critical temperature where the chiral condensate becomes zero in the real 
phase. Moreover, in the phase argL = vr of SU{2) gauge theory the chiral 
condensate is 1 for all temperatures 0. We shall show that this may only 
be valid in the GCE. Treating QCD within the CE leads to chiral symmetry 
restoration in all triality sectors. 

In the limit of vanishing quark mass the action for staggered fermions 
is chiral symmetric. Therefore, we turn now to the scheme with staggered 
fermions. For the effective action which includes effects of dynamical quarks 
in terms of mesonic field ^x'^x = c"x we use the effective action of Ref.[[l^ 
this action to the CE, we get 

2 



1 ^ roo 

exp[NVgf{L, (t)] = - E E / n d(rx exp[^L'' E SxSx+n + Vjm] 



k=l{s^},s^eZ{N) 

exp[-Nt E (^xV^-y^ay] [] [cosh NtE + ZkLs,] H 5 

x,y X X 

where (for d = 2) 

sinh E = m + ax 

and 

1 

^ n=l 



J2ax-Na 



.(14) 
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To get a relatively simple form for the EP we use the simplest mean-field 
approximation for mesonic fields Ox neglecting all fluctuations. Taking again 
the asymptotics for the two dimensional Ising model as described above we 
end up with the following expressions for the EP 

V^f (L, o) = ln(l -L'^)+\nH - Nta^ + lncosh(L/iJ) + 
2jL\tsinhL/Hy + -i'^L^[l + QiiaxihL/Hf - 7(tanhL/i7)^] + 0(7^) (15) 

for 7//^ ^ 1 and 

K^f (L, a) = ln(l -L^)+lnH- Nta^ + 2^L^ + F,{H) (16) 

for 7L^ ^ 1. The function Fc{H) is of the form 

F{H) = X/H + exp{-2X/H) + 2u^ exp{-AX/H) + 
u*[-exp{-AX/H) + Qexp{-QX/H) + exp(-8X/i7)] + O(m^), (17) 

where functions u and X were defined after Eg . ([T3|) . We have denoted here 

H = cosh NtE. 

These formulae are valid in the region X/ if -C 1 . In what follows we consider 
7 and Nt as independent variables and study EPs for different values of 7 and 
Nt. First of all, we observe that the EP in (|l^) for sufficiently small 7 has 
the maximum at L = 0. The EP in ([T6|) exhibits Z{2) symmetry breaking 
in the proper region of the convergence of the series in ( |17|) although, similar 
to the EP in ([T^), it is also center symmetric. Hence, with respect to Z[2) 
symmetry the behaviour of the system is close to the previous case with 
Wilson fermions. 

We suppose now that the EP (|l^) reflects all qualitative features of the 
theory in the region 7 < 1. This may be the case in the vicinity of L 
since the real expansion parameter is 7//^ and therefore the series in Eg . (|T5|) 
can still converge. To study the chiral symmetry behaviour we have plotted 
the EP (^) for several values of A''^ and 7 taking the quark mass m = 0. 
Fig. |l| shows the EP at 7 = 0.1 for Nt = 10 and for Nt = 4. Because these 
graphics are obviously symmetric with respect to Z{2) transformations we 
plotted them for positive L only. They demonstrate that the maximum of 
the EP is located at L = 0. One sees that at large Nt the chiral symmetry is 
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strongly broken (a 7^ 0). As Nt decreases (temperature increases) the chiral 
symmetry tends to restore and a ^ 0. Fig. ^ presents the EP (|l^) in the 
vicinity of the critical line for Nt = 2 and 7 = 0.8. One observes the main 
maximum at L 7^ and a = for these values of the parameters. 

In all cases for A^^ G [1 — 6] and 7 G [0.5 — 0.9] the spontaneous breaking of 
Z{2) symmetry is accompanied by chiral symmetry restoration since a = 
above the critical point. The phase transition is of first order. The important 
fact is that cr = in both sectors, i.e. for positive and negative values of the 
PL. Therefore, the chiral symmetry is restored in both sectors. This is in 
contradiction with the GCE description where only one phase exhibits such a 
restoration [0, §]. As 7 grows, the symmetric maxima of the EP are smoothly 
moving from each other to the values L = ±1 for the PL. In our approxi- 
mation it is impossible, however to establish reliably whether they can reach 
these values because in this region we are far away from the strong coupling 
where our EP can be trusted. We have shown this behaviour in Figs. |^-^ 
plotting contour graphics for Nt = 3. Fig. ^ presents the confinement re- 
gion with the only maximum at L = 0. Around 7 «i 0.6 two symmetric 
deconfining maxima with a = start to develop. Fig. ^ shows the critical 
point 7 = 0.68 where confining and deconfining maxima become degenerate. 
Fig. ^ shows the deconfinement phase where the maximum with L = is 
metastable. There exists a rather wide region in the (A^t, 7)-plane where we 
find the coexistence of the confining vacuum with cr 7^ and the two degen- 
erate deconfining vacua with a = 0. One of them is always metastable but 
physically acceptable because they have normal physical features in contrary 
to the so-called Z{N) metastable phase. This picture implies the existence of 
different " confining" and " deconfining" bubbles around the critical line with 
nonzero surface tension. One could interpret this vacuum structure as an 
indication of a mixed state in hot nuclear matter which has been discussed 
for more than ten years as a hyphotetical state in the phenomenology of the 
quark-gluon plasma but the existence of which, to our knowledge, has not 
been confirmed up to now from first principles. 

Let us briefly summarize our results and discuss some open questions. 

In this letter we gave a prescription for the calculation of the EP in the 
CE with respect to triality. We 1 that all our results are valid in = 3 + 1 
dimensions for the SU{3) gauge group though the analytical calculations 
in this case are more difficult. The EP in the CE demonstrates reacher 
phase structure in contrast to the corresponding potential in the GCE. Three 
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main points distinguish our approach from the conventional one: 1) The 
deconfinement phase transition manifests itself in the CE similarly to pure 
gluodynamics. 2) The chiral symmetry is restored together with spontaneous 
breaking of Z{Nc) symmetry and this restoration takes place in all n-ality 
sectors. 3) A reach vacuum structure around the critical line indicates the 
existence of a mixed state in hot nuclear matter. 

We showed these properties both for Wilson and for staggered fermions. 
The CE exhibits closer interrelation between deconfinement and chiral sym- 
metry restoration. This connection might be even tighter if we had not used 
the approximation L/H <^1. It could be very instructive in this respect to 
include a nonzero baryonic number in the present approach as well as a finite 
value for the quark mass. 



It has been argued in p, |T3[ that Z{Nc) global symmetry at finite temper- 



ature has a physical interpretation in terms of the charge of physical states. 
In our case this is the triality charge of fundamental quarks. The following 
question can be posed at this point: is the Z{2) symmetry really broken and 
L 7^ above the critical point or does the presented picture give only a for- 
mal way to discover a critical behaviour? The inconsistency may arise, since 
if the Z{Nc) symmetry is broken, there could be non-zero triality states con- 



tributing to the partition function [|1^. It would then intuitively seem that 
the projection onto zero triality states and the spontaneous breakdown of 
Z{Nc) symmetry are in contradiction with each other. The possible answer 
might be that the CE with respect to triality and projection of the physical 
observables onto zero triality states are different things: in the former case, if 
Z{Nc) is broken, there are nonzero triality contributions to the free energy at 
high temperatures, whereas in the later case such contributions are absent. 
However, we tend to the opinion that in the case of the CE description of 
finite temperature QCD non-zero triality contributions are all vanishing in 
the thermodynamical limit. That such a vanishing is not in contradiction 
with the Debye screening in the deconfinement phase has been discussed 
in An interesting approach is to consider the Hamiltonian formulation 
and to interpret all physics exclusively in terms of physical space variables 



11]. In this respect we reckon, that the application of the consideration of 



Ref. ||TT| to the theory with dynamical fermions described by the CE with 
respect to triality would give interesting results concerning critical behaviour 
of full QCD. We do not, however, think that Z{Nc) symmetry is unphysical 



as has been discussed in |]Tl| and |T^. It is really so, that Z{Nc) symmetry 
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acts as an identity on all physical states. The reason for this is Z{Nc) global 
transformations of zero temperature QCD together with the substitution 

for the fermionic fields \1/ determines the triality of quark states at finite 
temperature 0. In this sense, the Z{Nc) symmetry is a physical symmetry 
(i.e., if triality charges are physical charges). One has to be 0, nevertheless, 
treating this symmetry at finite temperature since it may occur that not 
only states having zero triality charge contribute to the partition function. 
The CE ensures that only states contribute to the partition function which 
transform identically under above transformations. 
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Figure 1: Effective potential in the confinement region for 7 = 0.1, iVt = 10 
above and Nt = A below. 
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Figure 2: Effective potential around the critical line for Nt = 2 and 7 = 0. 
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Figure 3: Contour Plot of the effective potential in the confinement region 
for = 3 and 7 = 0.5. 
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